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The most puzzling aspect of the glass transition observed in laboratory is an apparent decou- 
pling of dynamics from structure. In this paper we recount the implication of various theories of 
glass transition for the static correlation length in an attempt to reconcile the dynamic and static 
lengths associate with the glass problem. We argue that a more recent characterization of the static 
relaxation length based on the bond ordering scenario, as the typical length over which the energy 
fluctuations are correlated, is more consistent with, and indeed in perfect agreement with the typical 
linear size of the dynamically heterogeneous domains observed in deeply supercooled liquids. The 
correlated relaxation of bonds in terms of energy is therefore identified as the physical origin of the 
observed dynamic heterogeneity. 

PACS numbers: 64.70.Pf 



I. INTRODUCTION 

In the glassy phenomenology, it is the super- Arrehnius 
slowing down of the transport properties that is most 
striking 0, ■ We refer to the so-called fragile liquids 
0, which are distinguished by a highly temperature- 
dependent effective energy barrier, E c s(T), in their ther- 
mally activated expression for structural or a relaxation 
time 



T a {T) = ^ exp(£ off (T)/£; B T), 



(1) 



where Too ~ 10~ 13 s is a high-T relaxation time, and ks 
is the Boltzmann constant. The temperature variation 
of T a (T) for fragile supercooled liquids is described over 
a wide range of temperatures by the empirical Vogel- 
Fulcher (VF) equation 



T a(T) — Too exp[A/ (T — Tq 



(2) 



The apparent divergence temperature To appearing in 
Eq. J2J), is called the Vogel-Fulcher temperature, and of- 
ten found to be very close to the Kauzmann temperature 
Tk LJ where the configurational entropy of the liquid 
appears to vanish, if it were to stay in equilibrium [(|. 
Although Tk is not a rigorous lower-bound, it is a fairly 
good indication as to the lowest temperature a liquid can 
be supercooled Q- We note that Eq. |[2J automatically 
predicts a phase transition at T H It is the case 
that in finite dimensions, and short-range interactions, 
a diverging time is normally accompanied by a diverg- 
ing length. In fact, the large activation energy barrier 
E e fi(T g ) w 50 ksTg observed for weakly bonded frag- 
ile liquids at the laboratory glass temperature T g , is re- 
garded as an indication for the cooperative nature of the 
relaxation dynamics |ldj| . Several equilibrium theories of 
the structural glass transition, invoke an increasing static 
correlation length that diverges at Tk (~ Tq) [Tllll^ . ll3| . 
As the dynamics is activated, even a small correlation 
length leads to macroscopic values for r a exceeding the 
observation time: thus, the falling out of equilibrium of 



the liquid at the laboratory glass temperature T g (> To), 
which is conveniently defined as that temperature where 

T a (Tg) ~ 10 3 S i. 

The most puzzling aspect of this super- Arrehnius dy- 
namic slow down is that it is not seen to be accom- 
panied by any obvious change in the static structure 
of the liquid. The static structure factor (i.e. Fourier 
transform of the pair correlation function) obtained by 
X-ray and neutron scattering experiments reveals only 
marginal changes in the local density profile of polymers 
and molecular supercooled liquids [Tj|- In other words, 
an 'extended' short-range order of less than a nanometer 
in size begins to appear as the liquid is cooled toward T g . 
This occurs while the a relaxation time as obtained from 
viscosity and a.c. dielectric susceptibility measurements 
increases by as much as 12 orders of magnitude. It is 
reconciling the above two aspects of the structural glass 
transition, i.e. structure versus dynamics, which remains 
a major challenge in the condensed matter physics [l5|. 
Clearly, there is no detectable growing length associated 
with the density fluctuation. However, there are gen- 
eral arguments in favor of the fourth-order potential en- 
ergy density correlation function and associated thermal 
susceptibility, the structural specific heat, as the ideal 
quantities for studying the static properties in fragile 
(l6| and model glass- forming liquids |l7|. Although in 
experiments the contribution of the potential energy is 
difficult to separate from the kinetic contributions to the 
specific heat, it is accessible in computer simulations of 
model glass-formers such as prototype Binary Mixture of 
Lennard-Joncs (BMLJ) particles [H, d 

From another perspective, and unlike the simple liq- 
uids treated as homogeneous [2(j , deeply supercooled liq- 
uids are also distinguished by dynamic heterogeneous do- 
mains [2lll2^ |. a few nanometers across, with widely dif- 
ferent relaxation times (varying by five orders of magni- 
tude). Recent multi-dimensional Nuclear Magnetic Res- 
onance (NMR) measurements are found to favor a het- 
erogeneity siz e, £ het = 2-3 nm, i.e. 5 or more atomic 
diameters I2U 1231 . Results obtained from the fluctuation 
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theor y |2 4 | using the Heat Capacity Spectroscopy (HCS) 
data [25| are similar — though some times higher |2(j . As 
for their dependence on the temperature, more recent 
experimental procedures have discovered a growing dy- 
namic length accompanying the glass formation in col- 
loidal and molecular liquids [2]j. However, the origin of 
this dynamic heterogeneity remains unclear to date. Do 
dynamic heterogeneous domains correspond to any static 
correlation in the structure? The heterogeneity of time 
suggests possible heterogeneity in the structure. Thus, a 
knowledge of the typical size of the static and/or dynamic 
lengths, their temperature dependencies and correlation 
(if any) c an g o a long way in resolving the structural glass 
problem jlj 

In a previous work [lil |28| , we proposed a Bond Or- 
dering (BO) scenario for the glass transition in which the 
cooperative relaxation of bonds in terms of energy, uncor- 
related with density ordering or crystallization [29l l30l | , 
has been discussed. In this bond ordering picture, the 
structural component of the specific heat arising from 
the potential energy fluctuation, plays the central role 
as the thermal susceptibility associated with the static 
length for cooperative relaxation [l^. This static relax- 
ation length, £bo, is defined as that length-scale over 
which energy fluctuations are correlated. It may also be 
regarded as an operational definition for the average lin- 
ear size of Cooperatively Rearranging Regions (CRRs) 
discussed in the context of Adam-Gibbs-DiMarzio the- 
ory of cooperative relaxation [TH ^| . The scenario pre- 
dicts a growing (and possibly diverging) structural spe- 
cific heat with the lowering temperature for the fragile 
liquids, which has been in part corroborated by more re- 
cent Monte Carlo (MC) simulations of a Lennard-Jones 
binary mixture by Fernandez etal. They found 

that local potential energy fluctuations are in fact cor- 
related over distances much larger than the short range 
of the inter-atomic interactions. This is made evident by 
studying finite-size effects at normally inaccessible tem- 
peratures T < Tmct, where Tmo t is the apparent mode 
coupling transition temperature |3l|. Temperatures as 
low as 0.895 Tmct are accessed and studied in equilib- 
rium using a local swap MC dynamics that is insensitive 
to the 'cage-effect', and thus drastically reduces the equi- 
libration time. 

Our aim in this paper is to make evident the correla- 
tion between the static relaxation length, £bo> defined 
as the average length over which energy fluctuation is 
correlated, and the typical linear size of the dynamically 
heterogeneous domains, £ w, which close to T g is five 
or more atomic diameters |2ll |27| . This is done in an 
attempt to clarify the physical origin of dynamic hetero- 
geneity. 

The rest of this paper is organized as follows: In Sec- 
tion [HJ we briefly review the implication for the static 
correlation length, of the theories of structural glass that 
invoke an underlying ideal glass transition. Prominent 
examples are the Adam-Gibbs-DiMarzio entropy model, 
and discontinuous mean-field spin glass models that have 



a random first-order transition in the mean-field formal- 
ism. We also discuss the potential-energy-landscape view 
of dynamics in supercooled liquids in section [H] This 
paves the way in section IIIII for a phenomenological de- 
termination of the structural contribution to the specific 
heat, C s , and the corresponding static length of relax- 
ation, £bO; for the fragile supercooled liquids. The ther- 
mal behavior of £bo is derived from that of C s via an 
energy version of the Fisher scaling law. In section llVl 
we discuss the implication of this novel approach with 
regard to the physical origin of dynamic heterogeneity, 
and its apparent relation to the mean field models of 
structural glass transition. A brief summary of the main 
results appears in section Ivl 

II. IMPLICATION OF PREVIOUS 
THEORETICAL WORK FOR STATIC 
CORRELATION LENGTH 

A. Adam-Gibbs entropy model 

By considering the liquid to be made up of identical 
clusters of size z, each capable of undergoingindepen- 
dent rearrangement, Adam and Gibbs (AG) [Tj arrived 
at an activated expression for r a (T) where the effective 
energy barrier is given by E c ff(T) = z*(T)Afi. Afi is 
to be largely interpreted as the 'potential energy barrier' 
per atom against rearrangement in a cluster composed 
of z*(T) atoms, and taken to be a constant. z*{T) is 
the minimum 'critical' size of a CRR, allowing a tran- 
sition between two configurations. By definition, the 
configurational entropy of this minimal CRR is given 
by s* = fee In 2. It is easy to see that the ratio of the 
total number of atoms N to the size of minimal CRR 
z*, should be equal to the ratio of the total configura- 
tional entropy S c to the minimum entropy of a CRR s* : 
N/z*{T) = S c (T)/s* [13. Thus, z*(T) is expressed in 
terms of the configurational entropy of the system 

z*(T) = Nk B \n2/S c {T). (3) 

Furthermore, by assuming a hyperbolic form AC P — C /T 
for the excess specific heat of the liquid over the equilib- 
rium crystal [||, an approximate expression for the con- 
figurational entropy is obtained by the thermodynamic 
integration 

S C (T)= [ dTAC p /T = C{T-T K )/T. (4) 
Jt k 

On substitution into Eq. ©, this results in VF equation 
being recovered. We also note that z*(T) oc T/(T — 
Tk). Letting z*(T) ~ £agi & (minimal) static relaxation 
length is obtained 

Ug(T)=£ 00 [1-(T k /T)]- 1 / 3 , (5) 

where £oo is constant and should be regarded as the typ- 
ical inter- atomic spacing or an atomic diameter. 
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In fact the concept of the increasing size of cooperative 
regions with the lowering temperature, dates back to con- 
siderations by Jenckel in 1939 [22j. However, the ther- 
mal variation of cooperativity was formulated in terms 
of the atomic packing energy, while AG theory empha- 
sizes the role of configurational entropy. It should also be 
pointed out that the central results of the Adam-Gibbs 
theory could be obtained without recourse to the concept 
of configurational entropy |33j | with only the assumption 
of a diverging z*(T) as T — » T 2 , where T 2 is a second- 
order phase transition point. For instance it has been 
found that the ratio T g /T 2 — 1.28 ± 5.8% in agreement 
with the AG prediction T g /T K = 1.30±8.4%. As may be 
seen from Table I, the static length of cooperativity at T g 
obtained from Eq. (j^J in fact is too small to define coop- 
erative regions. Adam-Gibbs arguments give rearranging 
units with at most 10 atoms near T g |22|. Considerations 
based on mean-field models of structural glass transition, 
lead to bigger cooperative regions, as described next. 



B. Discontinuous mean-field spin glass models 

From the notion of droplet fluctuations together with 
the concept of random first-order transition observed 
in discontinuous mean-field spin glass models 0, l34| . 
a correlation length is predicted that varies as £mf ~ 
(T - T K )~ 2/d , where d is the space dimensions. To il- 
lustrate this mean-field result we follow the more direct 
approach presented in [34l |3f| . An statistical mechani- 
cal analysis based on replica method 0, Hill reveals in 
mean-field models of structural glass transition (p-spin 
models, Potts glass, etc.) that there are in fact two dis- 
tinct transitions at Td and Tr (Td > Tk), the first of 
which at Td is a dynamical transition akin to ideal mode 
coupling transition |3lj| . This dynamical transition at 
Td appears only in the mean-field, and is reduced to a 
crossover temperature in finite dimensions by the nucle- 
ation process or activated hopping between the multitude 
of the glassy metastable states, Af(T) m exp[Na(T)], in 
the range of temperatures Tr- < T < Td- Thus, the 
configurational entropy, S C (T) = fcslnA/"(T), maintains 
its extensive character for Tk < T < Td, as AT(T) is 
exponential in the number of particles N in this temper- 
ature range. However, as the thermodynamic transition 
point Tk is approached from above, the configurational 
entropy vanishes linearly, S c — Soo[(T/Tk) — 1)], and 
the system stays in a glassy configuration indefinitely. It 
should also be pointed out that the replica overlap or- 
der parameter changes discontinuously despite the lack 
of a latent heat. The replica approach |3J, |3(| further 
indicates that 



J\f(T < T K ) oc exp 
Af(T > T K ) oc exp 



T K 



T k -T 



A correlation length for the thermodynamic transition 
at Tk is obtained by matching Eqs. (JSJ in the critical 
region, giving N ~ (T — Tk)~ 2 , and letting TV ~ £^ F . 
Thus, 



Zmf(T)=£ 00 [1-(T k /T)} 



-2/d 



(7) 



-K 

V I — - 1 

T K 



(6) 



The correlation length exponent predicted is v — 2/3 in 
d = 3 dimensions. 

We note, however, that the cooperative regions ob- 
tained from Eq. Q involve at most 90 particles at T g . 
This is a significant improvement over the AG result, but 
falls short of the size of dynamic heterogeneous domains 
observed at T g , which typically comprise a few hundred 
atoms 26]. Another difficulty with this mean-field pic- 
ture is that its glassy behavior is greatly modified in the 
short-range versions when studied numerically in three 
dimensions [38l l39| . MC simulations of a short-range 
version of p-spin glass model in three dimensions |38l ] , in- 
dicate a continuous/second-order transition to the glassy 
phase accompanied by a diverging susceptibility or cor- 
relation length. The finite-size scaling results are found 
to be in favor of a correlation length exponent v = 1 (see, 
Fig. 3 in H^)- Similar behavior is observed in a short- 
range Frustrated Ising Lattice Gas (FILG) model [3^, 
which also is known to have a random first-order transi- 
tion in the mean-field. The correlation length exponent 
is found to be unity for FILG too. In view of Harris crite- 
rion for the relevance of disorder in a second-order phase 
transition [40l |. which requires v > 2/d, and the above 
observations, we believe that mean-field v = 2/3 should 
be regarded as a lower-bound estimate of the correlation 
length exponent for the fragile glass-forming liquids. In 
section IIIII we shall also argue for a correlation length 
exponent v = 1 for the fragile systems, using a different 
argument. 



C. Dynamics as an activated potential energy 
barrier crossing 

It is a long held view that in supercooled liquids, dy- 
namics is dominated by the topographic properties of 
the system's Potential Energy Landscape (PEL) 0, : 
long-time a relaxation is dictated by thermally activated 
crossing of the potential energy barriers separating differ- 
ent valleys of the potential energy surface $(ri, • • • , r^), 
which is defined over the 3N-dimensional configurational 
space of the liquid comprising N atoms. It is also ar- 
gued that activated transport over the potential energy 
barriers begins to dominate at low temperatures, where 
E e g(T) > 5ksT 0- This description of dynamics in 
terms of the (3N+l)-dimensional PEL, facilitates the 
study of collective phenomena in viscous liquids, and 
helps to unify in a simple way some of the static and 
kinetic phenomena associated with the glass transition 

E3- 

More recently, it has been demonstrated using molec- 
ular dynamics simulations of model glass-forming liq- 
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uids that the concept of activated hopping between 
whole super-structures of many PEL minima, called PEL 
Metabasins (MBs), is central to a quantitative descrip- 
tion of the long-time dynamics in glass-forming liquids 
[3 EEj . Here, the time evolution of the system is re- 
garded as a sequence of MB visits each with a residence 
time r. The mean residence/escape time from MB of en- 
ergy e is given by (r(e,T)) = T oc exp(E(e)/k B T), where 
e is defined as the energy of the lowest local minimum 
within the MB. In fact the activation barrier, E(e), is 
found to depend on the depth of MB in PEL, e, in a 
rather simple way E(e) re — e. The lower the 

e, the higher is the activation barrier E(e). A suitable 
average over the MBs visited by the representative point 
at a given temperature, gives the average residence time 
(r(T)), which is to be regarded as T a (T): 

W(t(T)> = J dep(e,T)exp{-E(e)/k B T) 

= exp(-£ off (T)/fc B T), (8) 

where p(e, T) de is the fraction of MBs within the range e 
to e+de visited by the representative point at T. Clearly, 
E c g(T) may be interpreted as some suitable average over 
the potential energy barriers E(e) (~ — e) encountered 
by the liquid at a given temperature ^5|. An spectac- 
ular demonstration of this assertion is the concurrence 
between the crossover to super- Arrehnius relaxation and 
the commencement of the variation with temperature of 
(e(T)> in an 80:20 BMLJ model liquid (see, Fig. 1 in 19]). 
In order to better illustrate the close correlation between 
E eB (T) and (e(T)}, in Fig. 1 we plot k B T ln^/r^) 
(= E c s(T)) against the average value of PEL minima, 
(e(T)), for the temperatures accessed using the data of 
Ref. [T3 • The correlation is rather impressive. This in- 
dicates among other things that the effective activation 
energy barrier embodied in the empirical VF equation, 
may also be taken as an estimate for the thermal varia- 
tion of (e(T)}, or that of the configurational energy den- 
sity ($(T))/AT, of the liquid [l6|: 



E e g(T) ~ -(<f>(T)}/N. 



(9) 



We use this result to estimate the temperature variation 
of the structural specific heat for the fragile glass-forming 
liquids. 



III. IMPLICATION OF BOND ORDERING 
SCENARIO FOR STATIC LENGTH OF 
RELAXATION 

In an attempt to clarify the physical origin of the dy- 
namic heterogeneity observed in deeply supercooled liq- 
uids, we expand on the concept of bond ordering [2^ and 
the associated correlation length ^]. As pointed out 
in section ITT1 £mf(T 9 ) is a significant improvement over 
£ag(T 9 ), but consistently smaller than the typical linear 
size of the dynamic heterogeneous domains, 6iot(T 9 ). As 



it turns out another static correlation length, ^bo(T 9 ), 
defined as the length-scale over which the potential en- 
ergy fluctuation is correlated, is in perfect agreement 
with the experimental observation. 

By bond ordering we refer to the correlated relaxation 
of bonds into their low-lying energy states, where inter- 
atomic bonds (as opposed to atoms) are treated as dis- 
tinct objects possessing internal degrees of freedom or 
energy levels. The length-scale over which the energy 
fluctuations are correlated thus defines the static length 
of relaxation, £bo- The fourth-order correlation function 
in terms of the local density p(r), which may also be in- 
terpreted as a two-point energy correlation function in 
terms of the local potential energy density (j>(r), is de- 
fined by g3 



G E (r) = (0(r)0(O))-(0) 2 
~ g(r 7 '6so) / 'r d ^' '. 



(10) 



Ge(t) is a fourth-order correlator as </>(r) oc p(r) 2 . The 
exponent 77', is related to the structural specific heat ex- 
ponent a, via a fluctuation-response equation 



1 



d d rG E (r) 



d d r/r 



d-2+1/ 



(11) 



where £bo is the static correlation length beyond which 
the correlation function rapidly vanishes. From Eq. Ijll|l . 
where the link between the structural specific heat C s 
and the static relaxation length £bo becomes apparent, 
an energy version of the Fisher scaling |48| is obtained: 



= (2-»7>. 



(12) 



Assuming 77' = for the fragile glass-forming liquids, 
which has been corroborated by the numerical studies of 
short-range versions of mean-field structural glass models 
IMHi, we get 



v = a/2. 



(13) 



The structural specific heat C' s is in fact the tempera- 
ture rate of change of the configurational energy. In view 
of Eq. 10 , it is approximated by 



dE cS 
dT 



(14) 



Using the effective energy barrier implied by equation 
©, i.e., E cS (T) = Ak B T/(T - T ), we have 



C s = Ak B T /{T-T ) 2 



(15) 



Eq. I|15fl implies a power law temperature variation for 
the structural specific heat of the fragile supercooled liq- 
uids, with an exponent a = 2. A critical power-law in- 
crease for the structural specific heat, has been rece nt!; 
reported in MC simulations of a BMLJ model liquid 
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The exponent v is thus given by v = a/2 = 1. Hence, 
with the effective potential energy barrier embodied in 
the empirical Vogel-Fulcher equation, we have 

£bo(T) = £ oo T/(T-T ). (16) 

In fact, the quantity F = T g /(T g —T ) (= £ BO (T s )/£oo) is 
some times regarded as a measure of the fragility of the 
liquid Q- Its range of values for intermediate to most 
fragile of the liquids is 3 < F < 7.5. Thus, more fragile 
the liquid, the larger is the static length of relaxation 
at T g . Using Eq. (Q2J, and T g /T = 1.28 ± 5.8%, the 
static length of relaxation at T g is determined to be in 
the range 3.4-7.0 inter-atomic spacings. This is indeed 
in perfect agreement with the typical linear size of the 
dynamic heterogeneous domains observed at T g [^H H3| . 

IV. DISCUSSION 

In Table I we compare the values obtained for static 
lengths £bo, £mf and £ag at T g . This is done us- 
ing the ratio T g /T = 1.28 ± 5.8% E^. We find 
£,Bo(T g ) = 3.4-7.0 atomic diameters in excellent agree- 
ment with £,het(T g ) > 5 atomic diameters. We also note 
that (energy) fluctuations of all linear sizes x < £bo 
must be possible, and indeed probable. (A similar, but 
stronger, effect arising from density fluctuations is be- 
lieved to be the cause of critical opalescence observed in 
fluids near their critical point |50j.) Thus, one expects 
heterogeneous domains with a distribution of relaxation 
times at T g , or not too far above To, such that the big- 
ger domains have longer life-times r. By assuming an 
activated form r ~ e x , and x to vary by a factor of 10, 
we have r values that vary by 4 decades. This therefore 
explains the existence of dynamic heterogeneities near 
T„ whose life-times differ by several orders of magnitude 
211. We see that the correlated relaxation of bonds in 
terms of energy, can be viewed as the physical origin 



of the dynamic heterogeneity observed in deeply super- 
cooled liquids. 

It also is a matter of considerable interest that the BO 
exponent, v — 1, is in perfect agreement with that ob- 
tained from three-dimensional simulations of models of 
structural glass that exhibit a random first-order transi- 
tion in the mean-field In fact, it appears that 
bond ordering scenario provides a physical interpretation 
for the replica symmetry breaking transition observed in 
theory and simulation, where members of coupled replica 
tend to sit in front of each other in the low tempera- 
ture glassy phase. Let us consider m coupled replica 
constrained to be in the same state. This resembles a 
liquid consisting of (super)molecular structures with ev- 
ery molecule composed of m number of coupled atoms 
[9j. This is an artifact of the attractive coupling among 
the replica, and can be viewed as a theoretical indication 
for the local ordering of bonds in the realistic systems. 
A detailed investigation of this apparent analogy will be 
presented elsewhere. 

V. SUMMARY 



To summarize, the exponent governing the thermal be- 
havior of the static relaxation length in fragile super- 
cooled liquids, v = 1, implies a stronger temperature 
dependence for cooperativity than previously thought. 
The static relaxation length estimated at T g is in per- 
fect agreement with the typical size of the dynamic het- 
erogeneous domains in deeply supercooled liquids. Fur- 
thermore, this larger static length also explains the wide 
distribution of relaxation times associated with the het- 
erogeneities, which is central to explaining the stretched 
exponential relaxation [5l| . and the decoupling of self- 
diffusion from viscosity in deeply supercooled liquids |52| . 
The correlated relaxation of bonds in terms of energy is 
therefore identified as the origin of dynamic heterogene- 
ity in fragile supercooled liquids. 
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FIG. 1: Dynamical activation energy E e g(T) is plotted 
against the average depth of PEL (e(T)) for temperatures 
accessed using the data of Ref. jlflj in order to illustrate their 
correlation. Solid line is the least squares fit to the data 
points. 



8 



TABLE I: The static correlation lengths £bo, £mf and £ag 
at T g axe tabulated for comparison. The observed linear size 
of dynamic heterogeneous domains |^],|27|, £,het{T g ) > 5, 
is in excellent agreement with the prediction based on bond 
ordering picture. 
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